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Abstract. A new molecular hydrogen cloud is found in the sub-damped Lya absorber [log Af(Hl) — 19.88 ± 0.05] 
at the redshift Zabs = 1.15 toward the bright quasar HE 0515-4414 (zom = 1.71). More than 30 absorption features 
in the Lyman band system of H2 are identified in the UV spectrum of this quasar obtained with the Space 
Telescope Imaging Spectrograph (STIS) aboard the Hubble Space Telescope. The H2-bearing cloud shows a total 
H2 column density iV(H2) = (S.Jtl'l) ^ ^O^*^ cm~^ and a fractional molecular abundance fn^ = (2.3t?;?) x 10"^ 
derived from the H2 lines arising from the J = — 5 rotational levels of the ground electronic vibrational state. 
The estimated rate of photodissociation at the cloud edge Jo < 1.8 x 10~* s~^ is much higher than the mean 
Galactic disk value, Imw ~ 5.5 x 10~^^ s~^. This may indicate an enhanced star-formation activity in the z — 1.15 
system as compared with molecular clouds at z ~ 3 where I ~ Imw- We also find a tentative evidence that the 
formation rate coefficient of H2 upon grain surfaces at z = 1.15 is a factor of 10 larger than a canonical Milky 



Way value, 71m w « 3 x 10 



The relative dust-to-gas ratio estimated from the [Cr/Zn] ratio is equal 



to = 0.89 ± 0.19 (in units of the mean Galactic disk value), which is in good agreement with a high molecular 
fraction in this system. The estimated line-of-sight size of L ~ 0.25 pc may imply that the H2 is confined within 
small and dense filaments embedded in a more rarefied gas giving rise to the z = 1.15 sub-damped Lya absorber. 

Key words. Cosmology: observations — Quasars; absorption lines — Quasars: individual: HE 0515-4414 



1. Introduction 

The most abundant interstellar molecule in the universe, 
H2, is currently observed not only in the Milky Way disk 
(e.g., Rachford et al. 2002) and halo (e.g., Richter et al. 
2003a), but also in the Magellanic Clouds (e.g., Tumlinson 
et al. 2002) and in more distant regions of the universe 
such as intervening damped Lya absorbers (DLAs) seen in 
spectra of background quasars (QSOs). The DLAs are the 
systems with neutral hydrogen column densities Ai"(Hi) 



They are believed to originate in pro- 



> 2 X 10^° cm 
togalactic disks (Wolfe et al. 1995). The systems with 
lower hydrogen column densities, 10^^ cm~^ < iV(Hi) 
< 2 X 10^° cm~^ , are formally called sub-DLAs. The sub- 
DLAs may also be related to intervening galaxies. At the 
moment there are known 9 molecular hydrogen systems 
detected in DLAs and sub-DLAs in the redshift range from 
Zabs= 1.96 to 3.39 (see Table 1). 



In this paper, we present results from the analysis of 
a new 10th H2 system detected at Zabs = 1.15 in the sub- 
DLA toward the bright quasar HE 0515-4414. This is the 
first detection of H2 with the STIS/HST at an interme- 
diate redshift - a cosmological epoch when an enhanced 
star formation rate (SFR) is observed in young galaxies. 
The SFR shows a peak at z ~ 1 over the redshift interval 
< z ^ 3 (see, e.g., Hippelein et al. 2003 and references 
therein) . 

Molecular hydrogen, being an important coolant for 
gravitational collapse of gas clouds at T ~ 10'^ K, is known 
to play a central role in star formation processes, and thus 
one may expect that the SFR and the fractional abun- 
dance of H2 are correlated. Studying H2-bearing cosmo- 
logical clouds leads to better understanding of the physical 
environments out of which first stellar populations were 
formed. 



Send offprint requests to: S. A. Levshakov 

* Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, 
which is operated by Aura, Inc., under NASA contract NAS 
5-2655; and on observations collected at the VLT/Kueyen tele- 
scope ESQ, Paranal, Chile, programme ID 066.A-0212 



2. Observations 

Spectral data of the quasar HE 0515-4414 (zom = 1-71, 
V — 14.9; Reimers et al. 1998) in the UV range were 
obtained with the HST/STIS (Reimers et al. 2001). The 
medium resolution NUV echelle mode (E230M) and a 
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Table 1. Molecular hydrogen abundances and dust contents k in cosmological H2-bearing clouds 



QSO 


Zahs 


log A''(Hi) 


logiV(H2) 


log /h2 


k* 


H2 detection 


0515 - 4414 


1.151 


19.88±0.05 


16.94+°:^? 




0.89 ±0.19 


this paper 


0551 - 366 


1.962 


20.50±0.08 


17.42l°;^g 


9 70+0.64 
^- ' °-0.90 


2.12 ±0.48 


Ledoux et al. 2002 


0013 - 004" 


1.973 


20.70±0.05 


19.84«:1° 




0.52 ±0.10 


Ge & Bechtold 1997 


1444 + 014 


2.087 


20.07±0.07 


18.30l°j^ 


1 /1O+0.38 


0.77 ±0.30 


Ledoux et al. 2003 


1232 + 082 


2.338 


20.90±0.10'' 


> 17.19'^ 


> -3.4r 


0.35 ±0.12 


Ge et al. 2001 


0841 + 129 


2.374 


20.95±0.10'= 


14.56^ 


< -5.98 


0.024 ± 0.005 


Petitjean et al. 2000 


0405 - 443 


2.595 


20.90±0.10 


18.161°:^^ 


^•^^-0.12 


0.17 ±0.06 


Ledoux et al. 2003 


0528 - 250 


2.811 


21.35±0.10'*''= 


'=18.22l°;23 


c 9 oq+0.25 
Z.OO_Q jg 


0.29 ±0.08 


Levshakov & Varshalovich 1985 


0347 - 3819 


3.025 


20.63±0.01 


14 61+° °^ 

iy:.01_Q Q2 


+0.02 
^- ' ^-0.02 


0.068 ±0.013 


Levshakov et al. 2002 


0000 - 2620 


3.390 


21.41±0.08-'^ 


n3.94t 


~ -7.2 


< 0.002 


Levshakov et al. 2000 



Note: Column densities A'^(Hl) and A'^(H2) are given in cm~^ ; ^tentative H2 identification; *photospheric solar 
abundances for Cr and Zn are taken from Grevesse & Sauval (1998), for Mg and Fe from Holweger (2001), 
see Eq.lOl. 

"High resolution UVES data reveal a few H2 subcomponents spread over ~ 700 km s"^ (Petitjean et al. 2002); 
''Srianand et al. 2000; '^Ledoux et al. 2003; '^Levshakov & Fohz 1988: '=M0ller & Warren 1993; 
^Prochaska & Wolfe 1999; "Levshakov et al. 2001. 
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Fig. 1. Typical example for a relative wavelength shift 
between successive orders. The H2 L2-0 P(3) line shows a 
difference of AA ~ 0.04 A (Az; ~ 5 km s~^ ). The lines at 
the bottom indicate the noise level. 

0.2" X 0.2" aperture provides a resolution power of A/AA ~ 
30, 000 (FWHM ~ 10 km s'^ ). The overall exposure time 
was 31,500 s. The spectrum covers the range between 
2279 A and 3080 A where the signal-to-noise ratio (S/N) 
per resolution element varies from S/N ~ 20 to ~ 5. The 
data reduction was performed by the HST pipeline com- 
pleted by an additional inter-order background correction 
and by coadding the separate sub-exposures. 

The spectral portion where the H2 lines occur suffers 
from a poor S/N ratio ( < 20). An additional problem 
arises from the limited wavelength accuracy. The MAMA 
detectors produce an absolute wavelength definition be- 



tween 0.5-1.0 pixel (2(7 limit as given by Brown et al. 
2002). For our data 1 pixel corresponds to 0.038 A. The 
spectral overlap of successive echelle orders allows to ex- 
amine the wavelength errors from order to order. Using 
well-defined line profiles we find relative wavelength shifts 
of 0.02-0.05 A (see Fig. 1 for an example). 

Additional echelle spectra of HE 0515-4414 were ob- 
tained during ten nights between October 7, 2000 and 
January 2, 2001 using the UV- Visual Echelle Spectrograph 
(UVES) installed at the VLT/Kueyen telescope. These ob- 
servations were carried out under good seeing conditions 
(0.47-0.70 arcsec) and a slit width of 0.8 arcsec giving 
the spectral resolution of A/AA ~ 55,000 (FWHM ~ 6 
km s-i ). The VLT/UVES data have a very high S/N ra- 
tio (~ 50-100 per resolution element) which allows us to 
detect weak absorption features. 

The high resolution VLT/UVES data reveal two nar- 
row components in the fine-structure C i lines associated 
with the sub-DLA at Zabs = 1-15 (de la Varga et al. 2000). 
The stronger component at Zabs = 1.15079 is separated 
from the weaker one by Au ~ 8 km s^^ , and shows 2.8 
times higher column density (Quast et al. 2002, hereafter 
QBR). Exactly at the redshift of Ci lines we identified 
more than 30 absorption features in the Lyman band sys- 
tem of molecular hydrogen H2. 

3. Measurements 

In this section we describe the measurements of the neu- 
tral hydrogen column density, metal and dust content and 
the H2 abundances in the Zabs = 1-15 sub-DLA. These 
values are well known to be physically related. The for- 
mation and maintenance of diffuse H2 in the Milky Way 
clouds is tightly correlated to the amount of interstellar 
dust grains, which provide the most efficient H2 forma- 
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Fig. 2. Continuum-normalized STIS spectra (histograms) of the quasar HE 0515-4414 (individual echelle orders 
marked by a, b, c, d, and e) and over-plotted synthetic H2 profiles (smooth lines) arising from the rotational levels 
J = to J = 5 of the lowest vibrational level v = in the ground electronic state X^T,'^. The noise level is shown by 
the dashed lines in each panel. The identified H2 lines are calculated for a two-component model with the components 
located at the measured redshifts of the Ci lines zi — 1.15079 and Z2 — 1.15085, having the widths 61 = 2.0 km s^^ , 
62 = 3.5 km s""'^ (Quast et al. 2002), and a column density ratio N2/N1 = 0.1. 



tion on their surfaces (see, e.g., Pirronello et al. 2000 and to be taken. Since the Doppler core of the Lya line is 

references therein). completely saturated, only the Lorentzian part gives in- 

formation about the line profile (Fig. (SJ. Moreover, the 
Lorentzian part is less pronounced than for typical DLAs 

3.1. Atomic hydrogen column density and metal and hence less distinguishable from the quasar contin- 

abundance uum. Therefore, we simultaneously optimised the contin- 

uum and fitted the spectral features using standard Voigt 

In order to estimate the column density of atomic hy- p^^gj^ ^^^-^^ technique. The continuum is modelled as a 

drogen contained in the sub-DLA, particular care has 
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Fig. 3. Part of the H i Lya forest showing the damped Lya 
system. The thin sohd and dashed hnes indicate the con- 
tinuum normahzed flux and the noise level, respectively. 
The thick solid curve represents the optimised Voigt pro- 
file model of the spectrum. The thick dashed curves show 
the individual profiles of the damped Ly a line correspond- 
ing to 70, 100, and 130 percent of the optimised column 
density N{Rl) = 7.6 x 10^^ cm'^ . The blue part of the 
damped line overlaps with some additional H i absorption 
components clearly identified by the presence of many as- 
sociated metal lines (Quast et al. 2003). The zero point of 
the radial velocity corresponds to the redshift z — 1.1508. 



linear function, and the Voigt function is calculated using 
the pseudo- Voigt approximation (Thompson ct al. 1987). 

Our optimized model (Fig. ^ reveals some additional 
absorption in the blue part of the damped Lya line at 
Aw ~ —420 km s^^ . This additional absorption is Hi 
Lya which is seen also in metal lines. The whole sub-DLA 
system is spread over 700 km s~^ (Quast et al. 2003). 
This line together with other narrow absorption lines seen 
in the wings of the damped Lya were included in the Voigt 
fitting. 

The derived column density of atomic hydrogen in the 
sub-DLA is log A^(Hi) = 19.88 ±0.05, where we estimated 
the standard deviation by varying the column density until 
the resulting model profile is apparently inconsistent with 
the observed data (Fig.|21). 

To measure the metal abundance in the main sub- 
component of the Znhs— 1.15 system we used the Znii 
AA2026,2063 hnes (Fig. Ej). The presence of dust grains 
in DLAs is usually estimated from the abundance ratio 
[Cr/Zn]^ assuming that Zn is undepleted (Pettini et al. 
1994). In our high S/N spectrum, only a weak Cm A2056 
line was detected at Av = km (Fig.^J). Other Cm 
lines (AA2062, 2066) are too weak to be visible. Their os- 
cillator strengths scale as /2056 : /2062 : /2066 = 1 : 0.74 : 
0.50 (Bergeson & Lawler 1993). 



^ [X/Y] = {X/Y)-{X/Y)q, where (X/Y) is the logarithmic 
value of the element ratio by number without reference to the 
solar value. Photospheric solar abundances {X/Y)q are taken 
from Grevesse & Sauval (1998) and from Holweger (2001). 
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Fig. 4. Parts of the UVES observations showing absorp- 
tion arising from the ions Cm and Znii (histograms). 
The solid and dashed lines represent our optimised model 
and its deconvolution, respectively. Mg i A2026 absorption 
from the same sub-DLA is indicated. The zero point of the 
radial velocity corresponds to the redshift z = 1.1508. 



The column densities for Znii and Cm were calcu- 
lated by Quast et al. (2003): logiV(Crii) = 12.01 ± 0.09 
and logA^(Znii) = 11.99 ± 0.02. We used these values to 
estimate the dust-to-gas ratio in Sect. 4.2. 

3.2. H2 column densities 

Molecular hydrogen at ^abs = 1.15079 is detected in the 
J = up to J = 5 rotational levels. At a spectral resolu- 
tion of ^ 10 km s~^ it is not possible to resolve the inter- 
nal structure in the H2 lines observed in the C i absorp- 
tion (see Fig. 1 in QBR). The H2-bearing gas may also be 
distributed over a wider velocity range as compared with 
C I which is easily ionised by UV photons in optically thin 
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Table 2. The H2 column densities for different rotational 
levels from the ^^abs = 

1.15 sub-DLA toward HE 0515-4414 



Level 


logiV(J) 


logiV(J) 


logiV(J) 




accepted 


min 


max 


J = 


16.47 


16.00 


16.70 


J= 1 


16.60 


16.00 


16.85 


J = 2 


15.85 


15.70 


15.95 


J = 3 


16.00 


15.90 


16.18 


J = 4 


15.00 


14.85 


15.30 


J = 5 


14.48 


14.30 


14.60 



zones. However, for a good approximation one can assume 
that H2 traces, in general, the volume distribution of C i 
since such correlation is indeed observed in the Milky Way 
(e.g., Federman et al. 1980). Therefore, in our H2 analysis 
we used a two-component model based on the observations 
of Ci by QBR. We note that the C i data were obtained 
with higher spectral resolution (FWHM ~ 5.5 km s^^ ) 
and considerably higher signal-to- noise ratio (up to S/N 
~ 130 for the parts of the spectrum with C l lines). 

Panels a e in Fig. 2 present cchellc orders of the STIS 
spectra of HE 0515-4414 (histograms) in the wavelength 
regions of the H2 Lyman 0-0 to 4-0 bands, together with a 
two-component Voigt profile fit of the data (smooth lines). 
It is seen that some of the identified H2 transitions are con- 
taminated by the Lya forest or blended with metals from 
different intervening systems. This hampers significantly 
the measurements of accurate equivalent widths and their 
analysis through the curve of growth. Moreover, the noise 
level (shown by the dashed line) is rather high for the 
available STIS data and this may explain why some of H2 
features are inconsistent with others. For instance, the ob- 
servational profile of L4P4^ in panel a differs from those 
of L3R4 (a), L3P4 (b), L2R4 and L2P4 (c), and L1R4 
(d). Relative strengths of the LORO and LORl lines from 
different cchclle orders (panels d and e) arc not consis- 
tent (LORO is partly blended with Fell A2383 at z = 0). 
The apparent depths of the close pairs L0P2 -|- L0R3 (e) 
and LOPl + L0R2 (e and d), as well as the single lines 
L1R2, L2R4 (c) are deeper than those calculated from the 
simultaneous fit to all H2 lines. 

Under these circumstances a standard x-square fit- 
ting cannot provide a statistically valuable measure of 
goodness-of-fit. To estimate model parameters we required 
that the calculated spectra were within 1 a uncertainty 
range for the majority of the unblended H2 profiles or 
their unblended portions which match the data. 

We tried to optimize a set of the H2 column densi- 
ties for the two-component model with the components 
located at the measured redshifts of the Ci lines zi = 
1.15079 and Z2 = 1.15085. The broadening 6-parameters 
were fixed at 61 = 2.0 km s"-'^ , 62 = 3.5 km s^^ (as de- 

^ Since all H2 transitions considered in the present paper 
arise from the ground electronic-vibrational state, we use a 
short notation like L4P4 which means L4-0 P(4) in the stan- 
dard form. 



duced from the Ci lines by QBR), but a column density 
ratio between the sub-components, N2/N1, was a free pa- 
rameter ranging from 0.03 to 0.36 (the latter corresponds 
to the Ci column density ratio found in QBR). 

For a given H2 component, the same 6-parameter was 
used independently of the rotational level. The column 
density in each J level was derived from several calcula- 
tions of the Voigt profiles with a fixed value of N2/N1 
and different N{J) which match the observational spec- 
tra. The limiting values of N(J) (an adjustable minimum 
and maximum) were chosen to estimate the uncertainty 
interval for column densities. 

All identified transitions from J = and 1 are optically 
thick, but the apparent central intensities of the LORO and 
LORl lines (d) are not zero (we consider the LORO line in 
panel e as corrupted by a bad merging of different spec- 
tra). These lines restrict A^(0) and ^''(l) by, respectively, 

5 X lO^s and 7 x 10^*^ cm^^ at N2/N1 = 0.03. On the other 
hand, we observe neutral carbon which is usually shielded 
in molecular clouds from the background ionising UV ra- 
diation by the H2 absorption arising from the J = and 
1 levels. An essential shielding in H2 lines occurs when 
A''(H2) > 10^® cm~^ . This gives us a hint at a possible 
range of 7V(J = 0, 1). 

For the lower value of N2/N1 = 0.03, the contribu- 
tion from the second H2 component is negligible, but the 
synthetic profiles are systematically narrower as compared 
with the data. The presence of the second component is, 
therefore, important. On the other hand, the maximum 
value of N2/N1 = 0.36 provides too wide synthetic pro- 
files even with N{0) = N{1) = 10^^ cm'^ . We found that 
with N2/N1 0.1 an optimal set of the H2 column densi- 
ties may be deduced. An example of such solution is shown 
in Fig. 2. The obtained results are given in Table 2. 

4. Discussion 

We investigate now the physical conditions in the /<^abs — 
1.15 H2-bearing cloud by considering the processes and 
parameters that balance the formation and dissociation 
of molecular hydrogen. Since our observations show a 
relatively high metallicity in this sub-DLA, [Zn/H] = 
-0.49 ± 0.10 [i.e., Z ~ (i-D^©], and the dust con- 
tent is approximately similar to the mean value for the 
cold gas in the Galactic disk {k ~ 0.9), we consider cat- 
alytic reactions on the surfaces of dust grains (HoUenbach 

6 Salpeter 1971) as the dominant H2 formation process, 
whereas ion- molecular gas phase reactions (Black 1978) 
are less efficient. 

The measured column densities can be used to esti- 
mate the kinetic temperature, Tkin, the gas density, nn, 
the photodissociation rate, /, and the rate of molecular 
formation on grains, R. However, in view of the large un- 
certainties in the column densities iV(0) and iV(l), we can 
only provide an order-of-magnitude estimate for these pa- 
rameters. 

Another obstacle in the H2 analysis is that the bal- 
ance equation is related to the space densities of H i and 
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Fig. 5. H2 rotational excitation in the Zabs = 1-15 sub- 
DLA toward HE 0515-4414. The logarithmic column den- 
sities N{J), divided by the statistical weight g{J) for each 
state J, are plotted against the excitation potential. The 
dashed lines represent fits from a theoretical Boltzmann 
distribution. The rotational ground states, J = and 1, 
fit on a line that is defined by a temperature of Tqi = 90 
K. The excited levels, J — 2 — 5, fit on a line with 
T25 — 400 ± 80 K, possibly indicating excitation mech- 
anisms through UV pumping and H2 formation pumping. 



H2. In case of homogeneous clouds one can assume that 
n(Hl)/n(H2) ^ A^(Hi)/iV(H2). However, this assumption 
may not be correct for DLAs where multiphase structures 
and complex profiles are usually observed. Observations 
show that with each step in increasing spectral resolution 
the profiles break up into subcomponents down to the new 
resolution limit. 

The sub-DLA at Zabs =1.15 reveals, for example, tran- 
sitions from neutrals and low ions (as C i, O i, C 11, Si 11) to 
highly ionized ions (as CiV, Siiv) spread over Av c± 700 
km s~^ (Quast et al. 2003) which implies that the neutral 
H2-bearing cloud (s) is embedded in a lower density, higher 
temperature gas. Neutral hydrogen H i can be spread over 
all gas phases that contain neutral gas with and with- 
out molecules. The H2 on the other hand may have a 
very inhomogeneous distribution in DLAs and concen- 
trate in small clumps (Hirashita et al. 2003). Thus, only 
a fraction of the total Hi may be relevant to the forma- 
tion of H2. In our Galaxy, for instance, 'tiny-scale atomic 
structures' (TSAS) and 'small-area molecular structures ' 
(SAMS) in the ISM are observed (e.g., Lauroesch & Meyer 
1999; Heithausen 2002). They show very high densities 
(nn ^ lO'^ — 10^ cm~^ ) and very small sizes (L ~ a few 
AU). 

To take this uncertainty into account, a scaling fac- 
tor (/< < 1 for the H i column density can be intro- 
duced. Following Richter et al. (2003b), who observed sim- 



n(Hi) 
n(H2) 



AT (Hi) 
A(H2) 



(1) 



where c stands for the cloudlet (s) where H2 is confined 
and 7V(H i) is the total neutral hydrogen column density 
along the sightline within the absorber. 

4.1. Kinetic temperature 

The kinetic temperature of the gas is usually estimated 
through the excitation temperature Tqi describing the rel- 
ative populations of the J = and J — 1 levels. This 
temperature is proportional to the negative inverse of the 
slope of the excitation diagram drawn through the points 
of the respective J levels in a plot \og[N{J)/ g{J)] versus 
E{J) shown in Fig. 5. Here, E{J) is the excitation energy 
of the rotational level J relative to J = 0, and g{J) is its 
statistical weight. 

Figure 5 shows that the value of Tqi is rather uncer- 
tain in our case because of large errors in N{0) and A^(l). 
Its mean value Tqi ~ 90 K corresponds to the excitation 
diagram shown by the dotted line, and its upper limit is 
about 270 K, which represents, probably, an upper limit 
for Tkin of the gas in the main sub-component of the Zabs — 
1.15 system. 

For levels with J = 2,3,4, and 5 the accuracy of the 
column densities is higher and we find T25 = 400 ± 80 K. 
The difference between (Toi)inax and (T25)inin is not sig- 
nificant and the points in Fig. 5 can be fitted, in principle, 
to a single excitation diagram. But the previous analysis 
of the fine-structure level populations of C i, where the 
most probable value for Tkin of 240 K was found (QBR), 
indicates that these two temperatures may not be equal. 
This is also in line with results on the H2 study in the 
Milky Way which revealed that single excitation diagrams 
fit usually only optically thin lines with 7V(H2) < 10^^ 
cm~^ (Spitzer & Cochran 1973). For higher column den- 
sities, there is 'bifurcation to two temperatures, depending 
on the J levels' (Jenkins & Peimbert 1997). 

4.2. Fractional molecularization and dust content 

According to our calculations presented in Sect. 3.1, the 
total H I column density in the main sub-component 
is A(Hi) ^ (7.6 ± 0.9) X 10^3 cni^^ . With A(H2) = 
(8.7^4 0) X 10^^ cm~^ , the ratio of H nuclei in molecules 
to the total H nuclei is 



/h2 — 



2iV(H2) 
A(H) 



= (2.3ti j) X 10 



(2) 



where Ar(H) = A^(Hi) + 2N{B.2)- 

Listed in Table 1 are the molecular hydrogen fractions 
in all known H2 systems. The values were derived in 
the standard way assuming the scaling factor (/) ^ I. This 
may imply that the listed molecular hydrogen fractional 
abundances are systematically underestimated. 



D. Reimers et al.: H2 at 2 = 1.15 toward HE 0515-4414 



7 



In Fig. 6 we compare these abundances with the dust- 
to-gas ratios, k, estimated from (Vladilo 1998): 

10[X/H]obs 



/x,ISM — /y,ISM 



10[X/Y]ob= 



1 



(3) 



where X and Y are two heavy elements with different de- 
pletions in dust. For all systems except Q 0405-443 and 
Q 1232-1-082 we used X = Znii and Y = Cm with their 
fractions in dust /zn,iSM = 0.587 ± 0.048 and /ci JSM = 
0.920 ±0.010 referring to the Galactic interstellar medium 
(Vladilo 2002, hereafter V02). Since column densities for 
Znii and Cm are not known for Q 1232-1-082, we used in 
this case X = Mg 11, Y = Fe 11 and /ms.ism = 0.715±0.056, 
/Fe,iSM = 0.939 ± 0.004 from V02, ahhough Mg and 
Fe may have not the same nucleosynthetic history. For 
Q 0405-443, the Cm abundance is not available and we 
used Fe 11 instead of Cr 11. The errors of the k values were 
calculated by applying error propagation method to the 
column density measurements quoted in the literature. 

Figure 6 demonstrates an apparent correlation be- 
tween and k in the range < fc < 1 which supports 
the assumption that molecular hydrogen abundances in 
quasar absorbers are governed by the dust content similar 
to that observed in the Galaxy. This conclusion rises the 
question: Why is the H2 detection in QSO absorbers in 
this case so rare (lower than 30% according to Ledoux et 
al. 2003) ? Following Hirashita et al. (2003), we suppose 
that a relative paucity of H2 observations in DLAs may 
be caused by a bias against finding H2 in dense molecu- 
lar clumps that have a small angular extent and thus a 
small volume filling factor. DLAs are mainly associated 
with diffuse clouds that have large volume filling factor 
and low molecular fractions, but they may also contain a 
small size dense filaments like the above mentioned TSAS 
or SAMS. Besides, low metallicity of the QSO absorbers 
can also significantly suppress H2 formation (Liszt 2002). 

One point in Fig. 6 (Q 0551-366) shows an unreal- 
istic high dust-to-gas ratio, about 2 times the Galactic 
value. This large value may be, probably, explained by 
systematic errors in the measurement of the Znii col- 
umn density. For instance, the relative abundance of Si, 
[Si/H] = -0.42 ± 0.11, differs significantly from [Zn/H] 
= -0.08 ± 0.12 according to Ledoux et al. (2002). The 
fraction of these elements in dust in the Milky Way 
is approximately identical, /si,isM = 0.691 ± 0.069 and 
/z„,isM = 0.587 ±0.048 (V02). 'At [Fe/H] = -0.90 ±0.11, 
Ledoux et al. measured [Si/Fe] ~ 0.5 which is in line with 
other observations (see Fig. 1 in V02). This means that 
[Zn/H] is most likely overestimated in the Zabs = 1.962 
system. 

4.3. H2 formation and photodissociation rates 

In equilibrium between formation on grains with rate co- 
efficient R (cm^ s~^) and photodissociation with rate / 
(s~"^), we may write that (Jura 1975b) 



O 
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Fig. 6. Relation between H2 fractional abundance 
plotted on a logarithmic scale and relative dust-to-gas 
ratio k (with respect to the mean Milky Way value) in 
DLAs and sub-DLAs. Indicated by numbers are the sys- 
tems found in spectra of QSOs 0000-2620 (1), 0841+129 
(2), 0347-3819 (3), 1232+082 (4), 0528-250 (5), 0405-443 
(6), 0515-4414 (7), 1444+014 (8), 0013-004 (9), and 0551- 
366 (10). A correlation between the two quantities in the 
range < fc < 1 is apparent (a large shift of fc for Q 0551- 
366 is probably caused by systematic errors in the data 
obtained by Ledoux et al. 2002, see text). 



where n — riu + 2nH2 1 ^^id /3o is the photoabsorption rate 
depending on the local UV radiation field (one may neglect 
the J dependence of (3 for an optically thick cloud since 
the photoabsorption rates from the levels J > 2 are low). 

To estimate the formation rate of H2 upon grain 
surfaces Rnn, we use approximation described by Jura 
(1975b). It assumes that (i) the levels J = 4 and J = 5 
are populated by direct formation pumping and by UV 
pumping from J = and J ~ I, (ii) the self-shielding 
in the levels J ~ and J = 1 is about the same, (iii) 
the upper levels J = 4 and J = 5 are depopulated by 
spontaneous emission (which is valid if njj < 10^ cm~'^ ). 
We do not consider additional rotational excitation of II2 
caused by a shock because restrictions on the gas density 
(riH ~ 100 cm~'^ ) and kinetic temperature (Tkin < 240 
K) set by the observations of Ci, Ci*, and Ci** (QBR) 
show that collisional excitation of the levels J — 4 and 
J = 5 is not significant. Using the cascade redistribution 
probabilities p4^o = 0.26 and P54 = 0.12, calculated by 
Jura (1975a), and assuming Tkin = 90 K, we can re-write 
Eqs. (3a) and (3b) from Jura (1975b) in the form 



/7iH2 — Rnn^i « 0.11/3o"-H2 7 



(4) 



[Rnn] 

and 

[Ruh] 



2.35 X 10 



9.15 X 10" 



iV(Hi)^ ' 



-9 ^^(5) 



iV(Hi) 



(5) 



(6) 
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By substituting numerical values in and © we obtain, 
respectively, [i?nH]c = {i-ltH) xlO-^'^(j)-^ s^^ and (3.6± 
1.2) X 1O~^'*0^^ s^^, which are consistent in the range 
[i?nH]c = (4.1 ± 1.7) X 10~^'*</)~^ s~^, and independent on 
the local value of (3o, as pointed out by Jura (1975b). 

To estimate the photodissociation rate at the cloud 
surface Iq, the shielding effect is to be taken into ac- 
count. The shielding factor, S, depends on line overlap, 
self-shielding of H2 , and continuum absorption. Lee et al. 
(1996) showed that these various factors can be well rep- 
resented by the H2 column density. They calculated the 
values of 5 as a function of A^(H2) for a turbulent veloc- 
ity of 3 km s~^ , which suits well for our case. From their 
Table 10 we find S = (1.928tJ;°^^) x lO^^ for N(H2) = 
(8.7^g y) X 10^^ cm^^ , respectively. This gives us a rough 
estimate of Iq = RniiN{}ii)/ SN(Il2) = 1-8 x 10"® s'^ 
(with the uncertainty of about 120%), or f3o ~ 1.6 x 10~^ 
s^^. The result obtained should be considered, however, as 
an upper limit on /3o since in our estimations we assumed 
that the H2 is one singe gas cloud. If, in reality, the H2 is 
inhomogeneously distributed among several cloudlets, the 
value of /3o should be lower. 

In the Milky Way, the mean value for /3o,haio ~ 
0.5/3o,disk = 2.5 X 10-1° s-i (e.g., Richter et al.' 2003b) 
and, thus, we may conclude that the H2 in the Zabs = 1-15 
sub-DLA is probably exposed to a radiation field with 
the intensity much higher than the mean Galactic value. 
For comparison, molecular clouds in the LMC and SMC 
also reveal 10-100 times more intense UV radiation field 
than the Galactic one (Browning et al. 2003). High ro- 
tational excitation of the H2 observed in the LMC/SMC 
gas and at ^abs = 1-15 is compared with Galactic data 
in Fig. 7. We note that the upper limit on the local UV 
field at 2:abs = 1-15 of about 100 times the Galactic value 
was independently determined from the analysis of the C i 
fine-structure lines by QBR. 

We can also estimate the photodissociation rate pro- 
duced by the intergalactic UV background (UVB) field on 
the surface of a cloud (Hirashita et al. 2003): 

/uvB 1.38 X IO-1V21 , (7) 

where J21 (in units 10^ erg cm^ "^s^i Hz^i str^i) is the 
UVB intensity at 1 Ryd averaged for all the solid angle. 
At Zabs= 1-15, J21 ~ 0.4 (Haardt & Madau 1996), and 
thus /uvB is of the order of 5 x 10-^'^ s-^, which is about 
4 orders of magnitude lower as compared with the derived 
value of Iq < 10^^ s-^. The presence of bright young 
stars in the Zabs = 1-15 sub-DLA is therfore required to 
maintain this high value of Iq. 

Star-formation activity is not, however, intense in the 
Il2-bearing clouds at higher redshift. For example, the UV 
radiation fields in the Zabs = 3.025 H2 absorber toward 
Q 0347-3819 and in the Galactic ISM are very much alike 
(Levshakov et al. 2002) . Other examples may be found in 
Ledoux et al. (2003). 
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Fig. 7. Total column density of II2 vs. excitation ratios 
N{A)/N{2) and A^(5)/iV(3), for LMC, SMC, Milky Way 
(data are taken from Browning et al. 2003), and sub-DLA 
at Zabs = 1-15 (dots with error bars). A systematically 
higher rotational excitation is observed in the Magellanic 
Clouds and in the sub-DLA as compared with Galactic 
data points. 

4.4. Gas density anc/ H2 formation rate coefficient 

We now consider constraints on the volumetric gas density 
riH and the II2 formation rate coefficient R stemming from 
the foregoing estimation of [i?nH]c ~ 4 x 10-^^0-"'^ s~^. 

A useful reference point is provided by the J = 2 level. 
The population of this level is more directly affected by 
coUisional processes, since it has a longer radiative lifetime 
as compared with J = 3 and other levels^. The critical 
density, , at which the probability of coUisional and ra- 
diative de-excitation of J = 2 are equal is 200 cm"'^ , if Tkin 
= 100 K, and = 80 cm'^ , if Tkin = 240 K (the col- 
lisional de-excitation rate coefficients qjj' are taken from 
Forrey et al. 1997). If nn > n^, coUisional de-excitation 
becomes important. 

To estimate nu, the grain formation rate of H2 should 
be known or vice versa, i? is a complex function of the gas 
and dust temperature and other poorly known parameters 
(HoUenbach & McKee 1979). According to their calcula- 
tions, i?MW ~ 3 x 10~i^ cm^ s'-^. Using this value, we find 
riH ^ 1300 cm^'^ , if < 1. However, for such large den- 
sity, the J = and J — 2 levels as well as the J — 1 and 
J — 3 levels should be in thermal equilibrium, which we 
do not observe. Moreover, the upper limit on the gas den- 
sity set from the excitation of Cl is 110 cm~^ (QBR). If 
we tentatively adopt njj ~ 100 cm^^ , then R ^ 10-Rmw- 

This value of R is larger than predicted in theoreti- 
cal calculations, but it is conceivable that R may vary in 
space since we observe considerable variations in the UV 
extinction among different clouds. A low grain formation 
rate of H2 (i? ^ 0.1i?Mw) was, for instance, recently esti- 

3 tj^2 = 3.4 X 10^" s, and tj=3 = 2.1 x 10^ s (Turner et al. 
1977) 
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mated in the LMC and SMC by Browning et al. (2003). 
Although our high value of R is consistent with the best 
determinations of upper limits to R toward 7 Peg, v Sco, 
A Sco etc. (Jura 1974), we cannot without further observa- 
tions conclude that this result is certain. To test whether 
or not the grain formation rate coefficient R exceeds the 
value of i?MWj a higher accuracy for the A^(0) and A^(l) 
column densities is needed to verify that Tqi ^ T25. 

For rtH ~ 100 cm^'' , the linear thickness of the H2- 
bearing cloud is small, L ^ 0.25 pc. Similar characteris- 
tics of molecular hydrogen small structures are found in 
intermediate-velocity clouds (IVC) in the Milky Way halo 
(Richter et al. 2003b). 

5. Conclusions 

We have analyzed the HST/STIS and VLT/UVES spectra 
of the quasar HE 0515-4414 and deduced the physical 
properties of the H2-bearing cloud embedded in the sub- 
DLA at Zabs — 1.15. The main conclusions are as follows: 

1. In the STIS spectrum of HE 0515-4414 we have identi- 
fied over 30 absorption features with the Lyman lines 
arising from the J = — 5 rotational levels of the 
ground electronic vibrational state of H2. These lines 
have exactly the same redshift as the fine-structure 
transitions in C i identified earlier by QBR in the 
UVES spectrum of this quasar. 

2. We find a total H2 column density of A^(H2) = 
(8.7140) ^ 10^® cm~^ and a molecular hydrogen frac- 
tion of = (2.3ti:?) X 10"^- 

3. From the measured ratios [Cr/H] — — 1.54± 0.11 and 
[Zn/H] = -0.49 ±0.10 we calculated the relative dust- 
to-gas ratio k — 0.89 ± 0.19 (in units of the mean 
Galactic disk value) in the molecular cloud (the ratio 
[Cr/Zn] is usually used to indicate the presence of dust 
in DLAs) . The derived H2 fractional abundance corre- 
lates with the dust content showing increasing fn^ with 
increasing k in cosmological molecular clouds. This in- 
dicates that catalytic reactions on the surfaces of dust 
grains is the dominant H2 formation process not only 
in the diffuse clouds in the Milky Way but in the DLA 
systems as well. 

4. Two excitation temperatures are required to describe 
the rotational excitation of the H2 gas: for the J = 
and J = 1 levels we find Tqi 90 K and the upper 
limit for the kinetic temperature of the gas of about 
270 K. For J = 2 - 5 we derive T25 = 400 ± 80 K. 

5. From the relative populations of H2 in the J = 4 and 5 
rotational levels we estimated the rate of photodissoci- 
ation at the cloud surface Iq < 1.8 x 10~* s^^ which is 
much higher than the mean Galactic disk value. Star- 
formation activity may be very intense in the close 
vicinity of the H2-bearing cloud, which is in line with 
the observed high SFR in galaxies at intermediate red- 
shift, z ~ 1. At higher redshift, z ~ 3, we do not ob- 
serve such intense UV fields in the molecular clouds. 



6. We also find that the formation rate coefficient of H2 
upon grain surfaces is probably 10 times higher as 
compared with the conventional value adopted for the 
Milky Way, Rmw « 3 x lO'^^ cm^ s'^. 

7. We find that in order to be consistent with the mea- 
surements of the population ratios of the fine-structure 
levels in Ci, the gas density in the molecular cloud 
should be ~ 100 cm~^ , which implies that the line-of- 
sight size of this cloud is small, L ^ 0.25 pc. Small sizes 
and a low level of the H2 detections in the DLAs favour 
calculations of Hirashita et al. (2003) who showed 
that the H2 may have very inhomogeneous distribution 
within these systems. The diffuse molecular hydrogen 
gas forms, probably, in small, dense filaments during 
the cooling and fragmentation phase in DLAs. 
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